Anomalous fluctuation regimes at the FFLO transition 
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Abstract 

Recently some experimental evidences have been obtained in favour of the existence of the inhomogeneous Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) superconducting state in heavy-fermion superconductor CeCoIns and organic superconductor 
A-(BETS)2FeCl4. However the unambiguous identification of FFLO state remains very difficult. We present the theoretical 
studies of the Gaussian fluctuations near the tricritical point (where the FFLO modulation appears) and demonstrate that 
the behavior of the fluctuational specific heat, paraconductivity and diamagnetism is qualitatively different from the usual 
superconducting transition. Special values of the critical exponent and the crossovers between different fluctuational regimes 
may provide a unique test for the FFLO state appearance. 



I. INTRODUCTION 

The long-standing hunt for the inhomogenous super- 
conducting state predicted by Fulde and Ferrell |l| and 
by Larkin and Ovchinnikov [2| (the so-called FFLO state) 
has recently known a revival in the study of fermionic 
cold atoms 0, 0L quantum chromodynamics color su- 
perconductivityl5| and heavy fermion superconductor 
CeCoIn 5 [f| 0, H|T The FFLO phase consists in a con- 
densate of finite momentum Cooper pairs in contrast to 
the zero momenta pairs of the usual BCS state. Hence 
the FFLO superconducting order parameter acquires a 
spatial variation. Such a state may be induced when 
a chemical potential difference is applied between two 
species of fermions with mutual pairing attraction. 

In superconductors the different chemical potentials 
are obtained by a Zeeman splitting h — (j,bH between 
spin up and down states, /is and H being respectively the 
Bohr magneton and an external magnetic field (in mag- 
netic superconductors h is the internal exchange field). A 
characteristic feature of the field-temperature phase dia- 
gram is the existence of a tricritical point (TCP) which 
is the meeting point of three transition lines separating 
the normal metal, the uniform superconductor and the 
FFLO state. For clean s-wave superconductors the TCP 
is located at T* = 0.56T c0 and h* = 1.07k B T c0 , T c0 be- 
ing the zero field critical temperature. Unfortunately, 
the FFLO state is only energetically favorable in a small 
part of the phase diagram, located at low temperatures 
T < T* and high fields, see Fig.??. Up to now, there is 
no compelling experimental evidence for this mysterious 
FFLO phase although some heavy fermion superconduc- 
tors like CeCoIns ml M, [H, organic superconductors like 
A-(BETS)2FeCi4 [9|, [nj, and rare-earth magnetic super- 
conductors like ErRh4B4 [ll[ are promising candidates. 

In the neighborhood of the TCP, the FFLO transition 
may be described within a modified Ginzburg-Landau 
approach (MGL). The MGL free energy functional, which 
follows from BCS theory with Zeeman pair breaking, 
was first obtained for s-wave superconductors [12|] and 
then gene ralized to d-wave pairing with/ without impu- 
rities [LI. The salient features of the MGL functional 
are the presence of higher order derivatives of the order 



parameter 'J than the usual g |V<J/| term, and the fact 
that the coefficient g of the latter term changes sign at 
the TCP being positive for h/2TTksT < 0.3 and negative 
otherwise. When the stiffness g becomes negative, the 
spatially varying configurations of the order parameter 
are favored leading to the FFLO transition. As a result 
of the anisotropy contained in higher order derivatives, 
the degeneracy over the FFLO modulation is removed 
even in the cubic lattice. Another source of anisotropy is 
related to momentum-dependent pairing. In the case of 
a ci-wave Pauli limited superconductor both lattice and 
pairing anisotropics are of primary importance for the 

fflo El [Tel. 

Beside these mean-field studies, the fluctuations have 
rarely been considered in the context of the FFLO tran- 
sition, although they may be accounted for in a quite 
simple and accurate manner within the framework of the 
MGL. Using the standard Ginzburg-Landau functional, 
the effect of fluctuations on the properties of the nor- 
mal state were extensively studied for the homogeneous 
BCS sup erconductors [13, EH and high-T c superconduc- 
tors [19( | . Though relatively small, thermodynamical and 
transport properties were actually measured in the Gaus- 
sian regime. 

Up to now, a detailed investigation of the fluctuations 
near the FFLO transition was missing. A first attempt 
to fill this flaw was realized recently by Marenko and 
Samokhin who addressed the issue of the quantum criti- 
cal fluctuations at the zero temperature FFLO transition 
pol ] . The aim of this paper is to perform the complemen- 
tary investigation of the fluctuations in the neighborhood 
of the tricritical point. 

In this Letter, we study the effect of the fluctuations 
close to the FFLO transition. We focus on the specific 
heat C and on the excess conductivity a, the so-called 
paraconductivity, induced by the superconducting fluctu- 
ations within the normal state. We have identified three 
regimes in the vicinity of the tricritical point (T*,h*) 
which differ by the nature of the soft fluctuation modes. 

Within the first regime, depicted as region I in Fig.??, 
the soft modes have the usual quadratic dispersion and 
are located around the origin of momentum space. This 
familiar situation yields the well-established power law 
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FIG. 1: Schematic (h,T) phase diagram around the tricrit- 
ical point (h*,T*). The solid lines represent the transitions 
between the normal metal, the spatially uniform BCS super- 
conductor, and the non uniform FFLO state. In region I 
where g > (resp. region II where g < ), the normal state 
behavior is dominated by uniform (resp. non uniform) super- 
conducting modes which propagate quadratically. In region 
III where g — > 0, the preponderant fluctuations propagate as 
e k ~ k 4 . Region III is defined by |T-T C | > g 2 /a-y. There 
are also crossover regimes (CR) separating the region III from 
I and II. 

behaviors (C,a) ~ (T - T c ) (d_4)/2 reviewed in [13, E, 

In contrast, the second regime represented as region II 
in Fig.??, is characterized by anomalous power laws in 
isotropic or weakly anisotropic materials. Here anoma- 
lous simply means distinct from (C, a) ~ (T — T c )^ d 4 ^ 2 . 
This discrepancy is related to the fact that the minima 
of the dispersion are attained on a finite dimensional sur- 
face, a sphere or a ring, see Fig.??.(c). Hence the avail- 
able phase space volume of the low energy fluctuations is 
larger than in the case of the standard BCS transition to- 
wards a uniform ground state, as it was first pointed out 
by Brazovski [21(. On the other hand, among the host of 
promising candidates for the FFLO state, crystal or gap 
symmetry effects are important. In such case the sphere 
or ring containing the low energy modes is partially re- 
moved and the low energy fluctuations are located within 
small islands centered at a few number of isolated points 
in momentum space, as it is shown in Fig.??. (a). Thus 
for large anisotropy, the anomaly disappears and one re- 
covers the usual power laws (C, a) ~ (T — T c y 4 ^ 2 . In 
summary the exponents in regime II depend drastically 
on the anisotropy of the system. 

The most interesting regime, region III of Fig.?? is 
at the immediate vicinity of the TCP where the FFLO 
and the uniform superconducting phases compete. The 
power laws of the specific heat C ~ (T — y c )( d ~ 8 )/ 4 an d 
paraconductivity a ~ (T — T c )^ d 6 ^ 4 are anomalous and 
quite universal since they pertain in both isotropic and 



anisotropic models. These predictions rely on the fact 
that the dispersion of the soft fluctuation modes becomes 
quartic within this regime. 

Indeed, these three regimes and the related crossovers 
between them may serve as a powerfull tool to identify 
FFLO state. 



II. GENERAL FORMALISM 

We investigate the fluctuation properties at the tran- 
sition between normal and FFLO states using a quit e 
general approach based on the MGL functional [lE, Il3|. 
Note that we assume a second order phase transition to 
the FFLO state. In the framework of isotropic model this 
transition is of second order for one and two dimensional 
systems while it is of first order in three dimensional case 
12]. However in the real superconductors the crystalline 
anisotropy (even a cubic one) is of primary importance 
and can modify this conclusion. Nevertheless, we may 
expect that in the case of weakly first order FFLO tran- 
sition, the fluctuation regimes studied here would be still 
observable. Such description is adequate near the tricrit- 
ical point where the most striking behavior is expected 
and where the wave-vector of the FFLO modulation is 
small. Looking for Gaussian fluctuations in the normal 
state it is sufficient to consider MGL functional only in 
the quadratic \^k\ 2 approximation 

d 

H i^ = J2 E ( a + 9 k f + njkftf) i*xi 2 (!) 

k i,j=l 

where V^k are the Fourier components of the complex 
order parameter and d the dimensionality of the sam- 
ple. The coefficient a — a(T — T c ), the coefficients g 
and jij depend on the magnetic field and temperature, 
T c (h) being the critical temperature for the second order 
transition between normal and uniform superconducting 
states. Although qualitative conclusions are quite gen- 
eral, explicit calculations are performed within the cu- 
bic anisotropy model defined by Yy = jSij 4- 777 (1 — dy), 
where 7 > 0. The case rj = 1 corresponds to the isotropic 
model since then the general fluctuation spectrum — 
a + gk 2 + Y.ijjijkfk' 2 depends only on k 2 = Sjfc^ and on 
/c 4 . Hence \j] — 1| is a measure of the deviation from the 
ideally isotropic material with a spherical Fermi surface. 

Integrating out the Gaussian modes 'i'k, one obtains 
the fluctuation free energy density 

L a 7t/cbJ 

k 

where L is the size of the system. The specific heat may 
be readily deduced from the free energy ([2]) via 

C=-T c %-a 2 k B T 2 c ^ y ■ (3) 

k k 
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may be rewritten as 



|J7-1|<1 



'7 = 1 



FIG. 2: Crossover between isotropic (c) and anisotropic (a) 
and (6) models in the regime \T - T c ] < ff 2 / a 7- Shaded areas 
represent the low energy modes leading to the various power 
laws. 



In evaluating the paraconductivity within the Kubo for- 
malism, one should use the time-dependent Ginzburg- 
Landau equation to obtain the current-current correla- 
tor at different times [l8[ . Moreover attention should be 
paid to the fact that, owing to the presence of high or- 
der derivatives jijk 2 k 2 in Eq.([2]), the expression for the 
current operator differs from the usual one. However the 
general expression fl9j relating the paraconductivity to 
the fluctuation spectrum, 



ne 2 ak B T c 1 ^ v i x 



4fi 



k 



■f. (4) 



still holds for the FFLO case although expressions for 
the velocities v^ x — dsk/dk x and for the current are very 
different from the usual quantum mechanical formula. 

Owing to the denominators in Eqs. (l3l4l) the main con- 
tributions to the specific heat and conductivity originate 
from the modes with the lowest energies £k- The above 
expressions lead to singular power laws which diverge at 
the transition with characteristic exponents. The expo- 
nents depend on the nature of the low energy modes: 
their location in momentum space and their dispersion. 
In the usual case of a uniform superconducting state, 
namely for g > and jij = 0, the low energy modes 
are concentrated around the center of the Brillouin zone 
and propagate with a quadratic dispersion. The result- 
ing power-law divergencies of the thermodynamical and 
transport properties are well-known [Tt], H, Eil ■ In the 
vicinity of the tricritical point, the situation is much 
more interesting since three different kinds of spectra £k, 
namely ([6|), |9]) or (fT3|) . lead to different power laws. 



III. ISOTROPIC MODEL AND LARGE STIFF- 
NESS 



We first consider the isotropic case r\ = 1. At low fields 
and high temperatures h/2nkBT < 0.3, the positive stiff- 
ness g leads to the standard spatially uniform supercon- 
ducting state (referred as BCS state on Fig.??). Con- 
versely when h/2irksT > 0.3, the stiffness becomes neg- 
ative and the fluctuation spectrum = a + gk 2 + -fk 4 



(•5) 



where q 2 = —g/2j and r — a — g 2 /Aj. Below T c = 
T c + g 2 /4: / ya > T c the normal state is instable towards 
the formation of the FFLO state with modulation wave 
vector q . Owing to the isotropy of the model, all direc- 
tions of q are degenerate. Above T c the minima of e£f ° 
are attained on a finite dimensional surface: a sphere for 
d = 3 or a ring for d = 2, whose radius is exactly q^. 
The one dimensional case is special in the sense that the 
minima consist in two points with coordinates ±go- 

At the vicinity of the sphere/ring defined by k = qo 
the approximation 



a (T - T c ) + A ig 2 [k - q y 



(6) 



may be used to evaluate Eqs. (|3l4p leading to the following 
singular behaviors for the specific heat 



fa\ 1/2 /lal\ (d ~ 2)/2 
C iso =A d k B T?(±) (M) (T-T c ) 

and for the paraconductivity 



■nA d e 2 k B T c 
M 



(2) 



1/2 



2 1 ) 



d/2 



(T - T c ) 



-■ill 



(7) 



(8) 



with Ad = l/4,l/8,l/87r for d = 1,2,3 respectively. 
Both d so and <7i SO diverge at the FFLO critical tem- 
perature T c with the same d-independent exponent in 
contrast to the usual uniform superconductivity case. 
In particular for the quasi-two and three dimensional 
cases, the divergencies of the fluctuation contributions 
near the FFLO transition are much stronger than the 
frp _ jy(d-4)/2 p 0wer j aw characterizing the usual su- 
perconducting transition. In this sense, the anomalous 
behavior of the Gaussian fluctuations is a signature of 
the FFLO state in such ideal isotropic systems. 



IV. ANISOTROPIC 
STIFFNESS 



MODELS AND LARGE 



The previous results (|71 8p rely deeply on the perfect 
isotropy of the model. Nevertheless real superconductors 
are anisotropic owing to either crystalline effect or <i-wave 
pairing. Hence the minima of the dispersion occur at a 
few number N a of isolated points which are located 
at a finite distance from the origin. These points belong 
to high symmetry lines of the crystal. More specifically 
in the model jij — "/Sij + 777 (1 — <5y) the wave vectors 

qQ Q ' are along the diagonals when 77 < 1 and along the 
crystalline axis for 77 > 1, whereas 7/ = 1 corresponds 
to the isotropic case. The relevant low energy fluctua- 
tions form small islands around the points , as shown 



3 



in Fig.??. (a). In the case of a very weak anisotropy 
\r) — 1| <C 1, these pockets may collapse within a finite 
dimensional surface and one recovers the isotropic model 
discussed above, see Fig.??. (6, c). 

Here we consider the opposite case of significant 
anisotropics \t) — 1| ~ 1 which is expected in real 
compounds. Then the favored modes belong to 
N a disconnected islands. Hence one may focus on each of 
the q$ obtained by solving 8,*EkO?o ) = and expand 
the fluctuation spectrum of (JTJ around q^ as 
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(a) 



(9) 



where A-"' = (1/2) (dieudjE^fal^') is the effective mass 

tensor for the fluctuating modes located around q^ . 
Owing to the symmetry of the crystal, the eigenvalues 
Ai of the tensors A,-°^ are the same for the N a islands. 

Evaluating ^ and (gj) with the spectrum © and col- 
lecting the contributions from the N a minima yield the 
specific heat 



(ah 



C„ 



N a k B T 2 c 



-d/2 



V / det Ai 



(T _ Tc) (d-4)/2 



(10) 



and the diagonal paraconductivity o~ a 
the plane xy 

N a e 2 k B T c ( (Ax + \ 2 )a d / 2 -' 



^det X 



(T - T c ) 



(d-4)/2 



(11) 

where Ai and A2 are the eigenvalues of \ff whose eigen- 
vectors span the xy plane. For the conductivity a zz 
one should replace Ai + A2 by 2A Z . In sharp contrast 
to the isotropic model results (|7I8|) . the temperature de- 
pendence of both Cani and o- ani becomes sensitive to the 
dimension with the same exponents as in the homoge- 
neous case [13, EH • Thus in the anisotropic model, the 
exponents themselves no longer provide a valuable hall- 
mark for the FFLO state. Naturally when r\ — > 1 one 
should recover the isotropic results (|7I8[) . The crossover 
from anisotropic to isotropic models is provided by the 
r\ dependence of det Ay . The latter determinant vanishes 
in the limit i] — > 1, as an example like det Ay = g 2 (l~n 2 ) 
in the quasi-two dimensional case. 



case) or in islands of size 5k centered at isolated points 
(anisotropic case). The width 8k, which is defined by 
the condition Aq^jSk 2 « r = a(T — T c ), must be smaller 
than go in order to have well defined shell or islands (like 
in Fig.??. (a)). Hence the common range of validity of 
Eqs. (l7l8|) and (|10|11|) is given by t <C 7<?o or equivalently 
\T — T c \ <C g 2 /aj- This corresponds to the area of the 
phase diagram which is sufficiently apart from the tri- 
critical point (large enough g) and close to critical line 
(small enough \T — T c \), schematically the region II on 
Fig.??. 

Conversely when 5k > q the relevant low energy 
modes merge towards the center of the Brillouin zone and 
the approximations ([6]) and ([9]) fail. Hence the specific 
heat and paraconductivity must be reevaluated in the 
regime \T — T c \ ^S> g 2 /aj. Then the physics is dominated 
by extremely soft fluctuation modes having a quartic dis- 
persion instead of the more usual quadratic one, yielding 
very special behaviors in the region III of Fig.??. 

Finally we investigate the crossover regime between re- 
gion I and III, see Fig.??. When g > the specific heat 
is given by 



knT, 



a Kb 



7 



d/4 



^n,(d-8)/i 



(l + x 4 + gx 2 1 ^forf)' 



(12) 

in the isotropic case, fid being the surface of the d- 
dimensional unit sphere. The parameter g controls the 
relative importance of quartic and quadratic terms. For 
large g the usual BCS results are recovered while the re- 
sult Eq. (TT4"|) is obtained in the opposite limit g — > 0. Both 
quartic and quadratic terms are equally important when 
g 2 r~j cry yielding to the same criterion \T — T c \ ~ g 2 /aj 
for the borderline between regions I and III than between 
regions II and III obtained above. 



VI. SMALL STIFFNESS FOR 
ISOTROPIC/ ANISOTROPIC CASES 



BOTH 



We now consider the model (TTJ in the regime 
|T — T c \ 3> g 2 /aj where the most salient features are 
expected. 

Approaching the tricritical point from the line defined 
by U/2-kT = 0.3, see Fig.??, one has g = and the 
spectrum of the fluctuations 



V. THREE FLUCTUATION REGIMES 

Let us check now the validity of the approximations ([6|) 
and ([9|) used so far in deriving the specific heat and para- 
conductivity. We aim to identify the region of the phase 
diagram where the previous results (|7I8|) and (|10llip ap- 
ply. The modes which participate to the singular parts 
of the specific heat and conductivity are spread within 
a shell of width 5k around the circle k = qq (isotropic 



4 = a(T-T*)+ ^ -,,/.f/,- 



(13) 



becomes purely quartic inducing the following anomalous 
temperature dependences of specific heat and paracon- 
ductivity: 



! \ 



d/4 



C* = k B T* 2 B d ( fx (77) (T - T *)( d - 8 )/ 4 (14) 
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with B d = 3^2/16,1/16, 1/16 V2tt for d= 1,2,3 respec- 
tively, and 



7r e 
4~fi 



k B T*C d 



(d-2)/4 



/ 2 fa) (T — T*) 



*\(d-6)/4 



( 15 ) 

with C* d = 3^2/32, 1/8tt, 5/96V2tt for ci = 1, 2, 3 respec- 
tively. The functions /i (77) and (v) are weakly depen- 
dent on the anisotropy factor and are equal to unity in 
the isotropic case rj = 1. Interestingly, the (/-dependent 
exponents differ from those which are typical of the uni- 
form superconductivity (region I on Fig.??) and of the 
FFLO transition at larger g (region II). 

Consequently, the regime \T — T c \ 3> <? 2 / a 7 (region 
III) and the related crossovers with regions I and II 
should reveal the existence of the tricritical point and 
FFLO phase. While they were derived along the line 
g = 0, the very unusual power laws C ~ (T - 7 1 *)( d - 8 )/ 4 
and a ~ (T — y*)(<i-6)/4 mus ^ emerge independently on 
the way the TCP is approached. 

Following the normal/FFLO transition line T c (h) to 
reach the TCP, one has g ~ r 1 / 2 where r = a{T* - T). 
We shall demonstrate that the expressions (|7I8I10I11|) 
applying in region II yield naturally the power laws 
(I14I15P accounting for the region III when T — > T* . 
First from Eqs. (|7ll0p . the specific heat is proportional 
to the product g( d ~ 2 )/ 2 -r~ 3 / 2 in the isotropic model and 
to g~ d / 2 T ( d ~ 4 )/ 2 in presence of sensible anisotropy, since 
the eigenvalues A.; of the effective mass tensor scale 
as g and detA ~ g d . Taking in mind that g ~ t 1 / 2 
along the transition line, one obtains the more diver- 
gent (T* — T)( d-8 )/ 4 behavior of the specific heat already 
predicted as Eq.(Q~J|) in region III for both isotropic and 
anisotropic models. We have used the fact that the FFLO 
and the second order uniform transition lines intersect 
themselves at the tricritical point. Furthermore simi- 
lar scaling on Eqs. (|8lllj) reveals that the conductivity in 
region II is proportional to g d / 2 r~ 3 / 2 for the isotropic 
model and to g( 2 ~ d )/ 2 -r( d ~ 4 )/ 2 for anisotropic models. In 
both situations, the substitution g ~ t 1 / 2 leads to the 
same power (T* — T)( d - 6 )/ 4 5 predicted above as Eq.lflg]). 
for the paraconductivity. In summary and more gener- 
ally, quartic terms in £k dominate the quadratic ones 
when \T-T C \ > g 2 /a^, and \T-T C \ < T c , see region 
III in Fig.l. 



VII. ORBITAL EFFECT 

Up to now the orbital effect was neglected to emphasize 
the influence of the Zeeman pair-breaking effect on the 
superconducting fluctuations. This is exact in ultracold 
fermionic atoms since there is no orbital effect at all for 
neutral objects. This is also a good approximation in su- 
perconductors with large Maki parameter which are the 
ones where the FFLO state is expected. Such a situation 
is generally encountered for i) layered superconductors 
under in plane magnetic field when the weakness of the 



interplane coupling quenches the orbital motion of the 
electrons (or for thin films of 3D superconductors), and 
ii) magnetic superconductors where the field originates 
from internal ordered magnetic moments. 

Beside these favorable situations, one should evaluate 
the effect of orbital pair-breaking on the results obtained 
previously, especially for three dimensional superconduc- 
tors under external magnetic field. Under strong mag- 
netic field, the electronic motion is described by the low- 
est Landau level wherein the kinetic energy is totally 
quenched except along the direction of the field. This 
provides a realization of an effective one-dimensional sys- 
tem. 

Coming from the BCS side with g > we evaluate how 
the fluctuation specific heat and paraconductivity scale 
as one approaches the tricritical point along the critical 
line T c (h), namely in the limit g — > + . For strong enough 
fields or sufficiently close to the transition line, the main 
contribution to the specific heat is given by the lowest 
Landau level as 



C = a 2 k R T 2 



2eB 



dk. 



9 k l 



(16) 



transforming the three-dimensional fluctuations into 
one dimensional fluctuations characterized by C = 
g~ l / 2 T c 3 / 2 ; f LUJc = 4geB/H being the Landau level spac- 
ing and t c = T+hLu c /2 = yielding the transition line un- 
der magnetic field. In a similar manner one may demon- 
strate that the magnetoconductivity along z-axis is given 
by a zz ~ g 1 ^ 2 Tc 3 ^ 2 . In sum one encounters a divergency 
of the specific heat while the magnetoconductivity a zz 
is not singular in the limit g — > + performed along the 
BCS transition line. 

Approaching the critical point along the line g = 
provides an interesting field-induced crossover. In the 
isotropic model the specific heat is then given by 



C = a 2 k B T 2 



2eB 



dk z 



2eB a 2 k B T 2 
h 7 V4 7 



► 7/4 



dx 

(l + a; 4 + Ca; 2 ) 2 



(17) 



(18) 



where r* = t 
terC-4( e B/fi)^7 



r )(2eB /H) 2 . The dimensionless parame- 
t* controls the relative importance 



of quartic an d qua dratic terms yielding a crossover field 
B* = (H/ejy/r*/^. A similar crossover between quartic 
and quadratic momentum dependence was already en- 
countered at the borderlines of the region III, see Eq. (fT2)) . 
Here the crossover is tuned by the magnetic field along 
the line g = 0. Similarly the conductivity 



e 2 7 
-irl — 



1/4 



5/4 



x dx 



(19) 



behaves as r* 5 / 4 for fields below B 
T *-3/2 f or i ar g er fields. 



(l + .x 4 + Ca; 2 ) 3 

fi(r*/ 7 ) 1 / 2 /e and 
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VIII. GINZBURG-LEVANYUK CRITERION 

Now let us address the issue of the validity of the Gaus- 
sian approximation. It is well known that Gaussian ap- 
proximation breaks down in the critical regime wherein 
the interactions between different fluctuation modes of 
the order parameter become sizeable. In the isotropic 
model, Brazovski demonstrated that due to the critical 
fluctuations the transition to the non-uniform state be- 
comes of the first order [2l[. Recently the same con- 
clusion has been obtained in the anisotropic model with 
isolated minima of £k [HJ • However performing a simple 
evaluation of the |^| 4 and \^f\ 6 interaction terms con- 
tained in the nonlinear MGL functional [13, EH, H3]> it 
can be proved that the relative width of the critical re- 
gion (in temperature or magnetic field) is given by the 
usual Ginzburg-Levanyuk parameter (T c /Ef) , Ef be- 
ing the Fermi energy [231 ] - More precisely and using our 
notations, the width tq of the FFLO critical regime is 
given by 



(T - T c )/T c = Tq 



Tc_ 

E F 



(20) 



when the fluctuations propagate quadratically, namely 
when \g/g \ > (T C /E F ) 4/3 (regions I and II of Fig.??). 
Here go is the value of the coefficient g far from the tricrit- 
ical point. Conversely when \g/go\ < (T c /Ef) 4 ^ 3 (region 
III of Fig.??) the critical region is somewhat broadened 
as 



iT-T, .'T, =- (! (-g- 



8/3 



(21) 



owing to the presence of quartic soft modes (region III 
of Fig.??). However even the largest value of the crit- 
ical region at g = is still very small, (T c /Ep) 8 ^ 3 ~ 
10~ 8 — 10~ 10 . Therefore the Gaussian approximation 



provides an excellent description of the fluctuation phe- 
nomena which may be observed at the superconduc ting 
FFLO transition. On the other hand the analysis [U, [22] 
may be relevant to the case of fcrmionic cold atoms where 
the ratio T c /Ep is not so small. 



IX. CONCLUSION 

We have demonstrated that the transition towards the 
FFLO state reveals very unusual and rich fluctuation 
regimes which may serve as a smoking gun of the FFLO 
state. Atomic Fermi gases provide a unique realization 
of an isotropic system while any real superconductor al- 
ways presents some degree of anisotropy owing either to 
crystalline effects and/or d-wave pairing. The fluctuation 
contribution vary non-monotonously near the tricritical 
point when we go from the uniform to FFLO state. In 
the vicinity of the tricritical point (region III) the very 
special exponents are expected for specific heat and para- 
conductivity divergencies owing to the presence of soft 
modes with unusual quartic dispersion. The presence 
of the tricritical point is signaled by very pronounced 
crossovers regimes which might be observed in experi- 
ments, by varying temperature or magnetic field. Finally, 
the fluctuation behavior of the magnetic susceptibility is 
also very peculiar at the FFLO transition and warrants 
detailed analysis [23| . 
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